On the Remarkable Spectrum of a Non-Hermitian Random Matrix Model 
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A non-Hermitian random matrix model proposed a few years ago has a remarkably intricate 
spectrum. Various attempts have been made to understand the spectrum, but even its dimension 
is not known. Using the Dyson-Schmidt equation, we show that the spectrum consists of a non- 
denumerable set of lines in the complex plane. Each line is the support of the spectrum of a periodic 
Hamiltonian, obtained by the infinite repetition of any finite sequence of the disorder variables. Our 
approach is based on the "theory of words." We make a complete study of all 4-letter words. The 
spectrum is complicated because our matrix contains everything that will ever be written in the 
history of the universe, including this particular paper. 

PACS numbers: ll.lO.Lm, ll.15.Pg, 02.10.Ud 



I. A CLASS OF NON-HERMITIAN RANDOM 
MATRIX MODELS 

Some time ago, Feinberg and one of us (in a paper 
to be referred to as FZ |^) proposed the study of the 
equation 



(1) 



where the real numbers Vj, are generated from some ran- 
dom distribution. Two particularly simple models were 
studied: (A) the r^-'s are equal to ±1 with equal prob- 
ability, and (B) = 6*^*= with the angle 9^ uniformly 
distributed between and 2w. 

Imposing the boundary condition ipo — ipN^-^ = on 
(|l|) we can write the set of equations as the eigenvalue 
equation 



(2) 



with the column eigenvector with components {ipi , ip2 , 
• • • , ipN, i^N+i} and Hn the {N + 1) by {N + 1) non- 
Hermitian random matrix 
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While quantum mechanics is of course Hermitian it is 
convenient to think of iJ as a Hamiltonian and as 
the non-Hermitian Schrodinger equation describing the 
propagation of a particle hopping on a 1-dimensional lat- 
tice. 

Some applications of non-Hermitian random Hamil- 
tonians include vortex line pinning in superconduc- 



tors y and growth models in population bi- 

ology^]. A genuine localization transition can oc- 
cur for random non-Hermitian Schrodinger Hamiltoni- 
ans §, 0, I §, 0, |ll|, |l|, |l| in one dimension. 

As mentioned in FZ, with the open chain boundary 
condition ipo = 'ipN+2 = the more general equation 



s/c+iV'fe+i + rk-it/Jk-i = Ei/jk 



(3) 



can always be reduced to (|I|) by an appropriate "gauge" 
transformation ipk — > Xkipk- Furthermore, applying the 
transformation ipk u~^iljk to (|l|) we see that if we 
change u^rk then the spectrum changes by — > 

uE. Thus, scaling the magnitude of the r^'s merely 
stretches the spectrum, and flipping the sign of all the 
Tfc's corresponds to rotating the spectrum by 90". 

It is also useful to formulate the problem in the transfer 
matrix formalism. Write (fil) as 



i^k+i 
tpk 



Ti 



fc-i 



i'k-i 



(4) 



where the transfer matrix Tk is deflned as the 2 by 2 
matrix 



Tk = 



E -Tk 

1 



(5) 



Define P = TnTn-i ■ ■ ■ T2T1. Then the bomidary condi- 
tion implies 



EPn + P12 = 



(6) 



The solution of this polynomial equation in E determines 
the spectrum. 

Since H is non-Hermitian the eigenvalues invade the 
complex plane. For model B, the spectrum has an ob- 
vious rotational symmetry and forms a disk (see Fig. 
which displays the support of the density of states). An 
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1 2 FIG. 2: Support of the density of states for model A for a 

" 4000 X 4000 matrix and a single realization of disorder. 



FIG. 1: Support of the density of states for model B. 



expansion of the density of eigenvalues around E ^ 
to very high orders in E has been given by Derrida et 
al. p4[ . This analytic expansion however cannot predict 
singularities in the density of states. 

In contrast, for model A the spectrum enjoys only a 
rectangular Z2 O Z2 symmetry. The first Z2 corresponds 
to E E* obtained by complex conjugating the eigen- 
value equation Hip = Eip. The second Z2 corresponds 
to E —t —E obtained by the bipartite transformation 
{--)''ipk- Remarkably, FZ found that the spec- 
trum has an enormously complicated fractal-like form. 
In Fig. H we plot the support of the density of eigenval- 
ues in the complex plane for a 4000 x 4000 matrix, for a 
specific realization of the disorder. In Fig. || we plot the 
support of the density of states in the complex plane for 
a 1000 X 1000 matrix, averaged over 100 realizations of 
the disorder. 

Contrasting Figs. ^ and ^ with Fig. |l|, one can see why 
it has been a challenge for mathematical physicists to 
understand the nature of the spectrum. 



II. BASIC FORMALISM 

In general, for random non-Hermitian matrices H, the 
density of eigenvalues can be obtained by 

1 d 

p{x,y) ^ --Q^G{z,z*) (7) 
where the Green's function is defined by 
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FIG. 3: Support of the density of states for model A for a 
1000 X 1000 matrix averaged over 100 realizations of disorder. 

with the bracket denoting averaging and z = x -\-iy (see, 
for example, ref. Q for a proof of these relations). Equa- 
tion (|^) follows from the identity 

d d 9 1 

— -logz=— - = 7r<5(x)%) (8) 
oz* oz oz* z 

where z — x -\- iy. Expanding G{z, z*) = 
i^^oZ-'=(^tri/'=) we see that {^tiH'') counts the 
number of paths of a particle returning to the origin in k 
steps. 

Evidently, for model A each link has to be traversed 4 
times. The spectrum of model A was studied by Cicuta 
et al. |l^ and by Gross and one of us (I6) by count- 
ing paths. In particular, Cicuta et al. gave an explicit 
expression for the number of paths. Recently, the more 
general problem of even- visiting random walks has been 
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studied extensively (see ref. jl^l)- In addition, exact an- 
alytic results for the Lyapunov exponent have been found 
by Sire and Krapivsky Q . 

In this paper we propose a different approach based on 
the theory of words. We will focus on model A although 
some of our results apply to the general class of models 
described by (^. 

One important issue is the dimensionality of the spec- 
trum of model A. In general, the spectrum of non- 
Hermitian random matrices when averaged over the ran- 
domness is 2-dimensional (for example, the spectrum of 
model B). Many of the authors who have looked at model 
A beheve that its spectrum, as shown in Figs. || and ||, 
is quasi-O-dimensional: the spectrum seems to consist of 
many accumulation points. Here we claim that the di- 
mension of the spectrum actually lies between 1 and 2 
dimensions, in the sense described below. 



III. DISTRIBUTION OF THE 
CHARACTERISTIC RATIO 

Consider the degree {N + 1) characteristic polynomial 
'-^N+iiz) = det(i?7v + zl^+i), where In+i denotes the 
(A^ + 1) X (A^ + 1) identity matrix. We easily obtain the 
recursion relation 

An+i{z) = zANiz) - rN^N-iiz) 

with Aq{z) = 1 and Ai(z) = z. Note that this second or- 
der recursion relation can be expressed in terms of trans- 
fer matrices as 



AAr+i(z) 

An{z) 



-rN 




Ajv(^) 
An-i{z) 



(9) 



very similar to those defined for the wave function ip, 
with the transfer matrix given by 



N 



-rN 




Following Dyson and Schmidt ^ 
characteristic ratio 

Ak{z) 
Afc_i(z) 

which satisfies the recursion relation 



(10) 



we consider the 



Vk 



(11) 



(12) 



with initial condition y\ = z. 

The hope is that, while the characteristic polynomials 
Afc(z) obviously changes dramatically as k varies, the 
characteristic ratio yk{z) might converge asymptotically. 

From the definition in eq. ( pH ) , it is clear that a point z 
is in the spectrum of the N x N matrix iff yN{z) = and 
yN+i{z) = oo. Thus a point z belongs to the spectrum 



if the corresponding set of variables {yjv} is unbounded, 
that is, if the probability of escape of the variable 2/„ to 
oo is finite. As is shown in jl^, this condition is sufficient 
to determine the spectrum along the real axis (z S R), 
but is insufhcient in the complex case. 

Let Pkijjk) be the probability distribution of yk (note 
that yk is complex). Then 



Pk+iiyk+i) 



d^yk Pk{yk)S^^Hyk+i -z + —) 

Vk 



where the brackets denote the average over the disorder 
variables {r^}. In the thermodynamic limit k oo, 
under fairly general conditions pl[ it can be shown that 
the probability distribution has a limit P{y), called the 
invariant distribution, which is determined by the self- 
consistent equation 



P{y) 



dhP{t)6^^\y-z+^-)^ 



-y\i \-y 



where we have used the fact that ^'■^•'(/(z)) near a zero, 
zo, of f{z) is given by 5^^\z - zq) 1^1"^ 



For model A, we obtain the amusing equation 



P{v) 



2 \ z- 



y 



z-y 



Pi: 



z-y 



(13) 



This type of equation has been studied extensively for 
the real case in [|2], Q . It can be shown that it can 
be solved by the Ansatz 



P{y) = Y,a,6(^\y- 



b,iz)), 



(14) 



where the 5j's depend on z whereas the aj's don't. Note 
that the index j is not necessarily an integer and can refer 
to a continuous set. In addition, it can be shown that 
the bj (z) are the stable fixed points of the product of any 
sequence of transfer matrices U (see the next section). 
Plugging (14) into (HSh we have 



j 3 

+5^^\y^ 



0, 



Since the right hand side has twice as many delta function 
spikes as the left hand side, for the two sides to match 
we expect that, in general, the index j would have to run 
over an infinite set. 

For a given complex number z, we demand that the 
two sets of complex numbers {bj} and {z— 1/bj, z + l/bj} 
be the same. This very stringent condition should then 
determine the {bj{z)}. 

To see how this works, focus on a specific 6i (Since 
the label j has not been specified this can represent any 
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bj). It is equal to either z — I/62 or z + I/62 for some 
62. But 62 must in its turn be equal to either z — I/63 or 
z + 1/63 for some 63. This process of identification must 
continue until we return to bi. Indeed, if the process of 
return to bi occurs in a finite number L of steps, then 
it will repeat indefinitely (since the system is back at its 
starting point 61). It is this infinite repetition which gives 
a finite weight to the (5— function at 61. By contrast, if 
the number of steps needed to return to the initial point 
bi is infinite, then the weight associated with this point 
vanishes, and it will not be present in the spectrum. 

We thus conclude that the support of the distribution 
of P{y) is the closure of the set of all the stable fixed 
points of the product of any sequence of transfer matri- 
ces U. We also conclude that the support of the density 
of states of the non-Hermitian matrix, in the thermody- 
namic limit, is given by the zeroes of any stable fixed 
point: bj{z) = 0. What is important to notice is that the 
Qj are independent of z and depend only on the length and 
of the word. Thus, we conclude that the set of complex 
numbers {bj{z)} is determined by the solution of an in- 
finite number of fixed point equations. 



IV. THE THEORY OF WORDS 

It is useful here to introduce the theory of words. 
A word w of length L is defined as the sequence 
{wi,W2, ■ ' ' iWl} where the letters Wj = ±1. In other 
words, we have a binary alphabet. Let us also define the 
repetition of a given word a specific number of times as 
a simple sentence. We can then string together simple 
sentences to form paragraphs. 

For a given word w of length L, consider a function 
fL{b;z,w) to be constructed iteratively. For notational 
simplicity we will suppress the dependence of /l on 6, z, 
and w, indicating only its dependence on the length L of 
the word w. The iteration begins with 



and continues with 



(15) 



We define f{b;z,w) = f^. 

The set of complex numbers {bj{z)} are then deter- 
mined as follows. Consider the set of all possible words. 
For each word w, determine the solution of the fixed point 
equation 

b = f{b;z,w). 

By considering small deviations from the solution, we see 
that the solution is a stable fixed point only if 



df{b;z,w) 



db 



< 1. 



(16) 



The set of all possible words w generates the set of com- 
plex numbers {bj{z)}. In other words, b is determined by 
a continued fraction equation, since 



f{b;z,w) 



We see that has the form 



ajb- 



wl 



(17) 



with the polynomials aj and (3j determined by the recur- 
sion relations 



za.; 



zA- 



(18) 



(19) 



with the initial condition ao = 1, ai = z, /3o — 0, and 
/3i = — wi. We notice that aj and Pj satisfy the same 
recursion relation as that satisfied by A, with the corre- 
spondence Wj+i ^ rj. Note also that ( [L8[ ) and (|l9|) can 
be packaged as the matrix equation 



1 



a,- 



, (20) 



where the transfer matrix Uj defined in the previous 
section appears. This is closely related to the trans- 
fer matrix formalism discussed earlier. Indeed, defin- 

aj-i /3j_i 

z —Wi 

1 

can also be characterized by a matrix 



ing Wj 
Wi = 



we have the initial condition 
Hence a given word w of length L 



W = 



a (3 
7 S 



(21) 



where for convenience we have written a = ul, P 
7 = ttL-i, and S = 

For a given word w, the fixed point value b is deter 
mined by the quadratic equation 



b = 



ab + P 
jb + 6' 



(22) 



which is the fixed point equation of the homographic 
mapping associated with the matrix M. The geo- 
metric interpretation is clear: The matrix W acts on 
2— component vectors w, and we ask for the set of v such 
that the ratio b of the first component to the second com- 
ponent is left invariant by the transformation. In other 
words, we look for the projective space left invariant by 
the transformation W: the fixed point value b defines the 
direction of the invariant ray. 
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Hence b is given by 



b ^ 



P{z) ± ^Q{z) 
2R{z) 



(23) 



with P,Q,R polynomials of degree 2L in z, where L de- 
notes the length of the word w. Explicitly, 



Pl(2) = aL- 



Ql{z) = (aL - Pl-i? + ^aL-iPi 



and 



Rl{z) = aL-i- 



(24) 



(25) 



(26) 



The stability condition ( p^ ) determines which root of 
( p3| ) is to be chosen. 

We will see shortly that b{z) determines the spectrum. 
Anticipating this, we see that if we form a compound 
word by stringing the word w together twice (for exam- 
ple, the Japanese word "nurunuru") then we expect the 
contribution to the spectrum to be the same. But given 
the preceding discussion, this is obvious, since if a ray is 
left invariant by W, it is manifestly left invariant by . 



V. DENSITY OF EIGENVALUES 

Once we have determined bj{z), how do we extract 
the density of eigenvalues? The eigenvalues {A^} of the 
matrix H are given by AAr(z) = Il^^{z — Xi). From ( [ll| ) 

we have J2k=i^ogyk ^ logAAr(z) = Eili log(z - Xi), 
and thus 



A.) 



(27) 



Using the identity (|8|) we can differentiate the right hand 
side of (|27| ) to obtain the density of eigenvalues in the 
complex plane 



- 1 9 d 
^ TT dz* dz 

Plugging in our solution 

P{y) 



E 



(PyP{y)logy. 



a,S^^\y~b,) 



rfr^y aAogb, 



we finally deduce that p — 

„ 1 9bj 

Since the aj do not depend on z, the spectrum is de- 
termined by the zeroes of the fixed point solutions bj{z). 

We see from (23) that the density of eigenvalues is 
given as a sum over j of terms like 



_d_ 

dz* 



P{z) ± ./Q{z) 



[P'(z)± 



2./Q{z) 
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FIG. 4: Support of the density of states of the periodic word 
{+ + -}• 



Thus the spectrum consi sts of isolated poles given by the 
zeroes of R and P(z) ± \jQ{z), and of the cuts of \jQ{z), 
and is made of isolated points plus curved line segments 
connecting the zeroes of Q{z). 

Contrary to what some authors have believed, 
the spectrum is not 0— dimensional, but (0 -I- 1 -I- 
(5)— dimensional, with 5 < 1. Each word gives rise to 
a line segment, and words which differ slightly from each 
other gives rise to line segments near each other. Indeed, 
given a word w, it is possible to construct a word w' with 
a spectrum as close to that of w as desired. For that pur- 
pose, we may construct w' as w' = n/^vw^'^ where v is 
any "corrupting" word, and the two lengths li and I2 are 
sufficiently long. Indeed, in terms of transfer matrices 
and invariant rays, we see that u/^ acting on any ini- 
tial ray brings it close to the stable invariant ray of w. 
Then the direction of this ray is corrupted by u, but it 
is brought back arbitrarily close to the invariant ray of 
w by applying the transfer matrix it;'^ , provided that li 
is large enough. Presumably (although this remains to 
be proven rigorously), the spectrum associated with the 
corrupted word w' can be made as close as we want to 
that of w. We have thus this property that for any word 
w, there is a word w' generating a spectrum as close as 
we want to that of w. In Figs. ^ and || we plot the eigen- 
states of a word w = {+ H — } and the spectrum of the 
word w' = + +}w^^. We see that the two spectra 

are very close. 

As is clear from this discussion, the spectrum is indeed 
"incredibly complicated." 



VI. WORDS AND SPECTRAL CURVES 



We content ourselves by focusing on the cuts of yQ{z). 
Since for a word w of length i, Q{z) is a polynomial of 
degree 2L with 2L roots, it gives rise to L curved line 
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FIG. 5: Support of the density of states of the periodic word 
{+ -\ — } corrupted. 

segments. The curves are given by the condition 

ImQ(z) = 0. (28) 

(The sign of Re Q{z) depends on the choice for the square 
root branch cut.) Given a word w of length L, the cor- 
responding spectrum must be invariant under a cychc 
permutation of the letters 'Wi,W2, ■ ■ • , wl, namely under 

Wi W2, W2 ^ W3, - ■■ ,Wl ^ Wl. 

As an example, for the 3— letter word w = {+ H — }, 
R{z) ^ z"^ + 1 and Q{z) ^ z^ - 2z^ + z"^ + A, which has 
roots aX z = ±i, z = ±(a/7-z)/2, and z = ±{\/l + i)/2. 

It is now clear what the words correspond to "physi- 
cally": a matrix H with r^'s given by an endless repeti- 
tion of (-1-1, 4-1, —1) has a spectrum given by a straight 
line connecting z = ±i, and two algebraic curves connect- 
ing z = {^/^ + ^)/2 to z = (\/7-i)/2 and z = -(\/7-hz)/2 
to z = — (-^7 — «)/2, plus poles at z = ±i. Notice that 
the two poles are buried under a cut. In Fig. ^we show 
the spectrum associated with the word {+^ — } together 
with the spectrum of a random 1000 x 1000 matrix. 

We now give a complete study of all 4~letter words. 
The polynomial is easily found to be 

Qi{z) = z* - 2sz'^ + (s^ + 2k)z'^ - 2sKz2 -I- uj"^ (29) 

with s = Wl + + W2 + Wi, K = W1W3 + W2W4, and 
u! = W1W3 — W2W4. The condition ( |2^ ) for the curves in 
the spectrum reduces to 

xy - + (2;* + ~ Qx^y^ + s{y^ - x^) + fc) = 0. 

(30) 

There are only three non-trivial 4— letter words, namely 
w = {+ + H — }, {+H }, and {H }. Their contri- 
bution to the spectrum of H together with the spectrum 
of a random 4000 x 4000 matrix is shown in Fig. |^. 

In Fig. ^ we show the contribution of all one, two, 
three, and four letter words to the density of states. 
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FIG. 6; Support of the density of states of the word {+ H — } 
over the spectrum of a random 1000 x 1000 matrix 
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FIG. 7: Support of the density of states of the periodic words 

{-|- + + — }, {+ H } and {H } over the support of the 

spectrum of a random 4000 x 4000 matrix. 



Thus, an by iV matrix H with r^'s given by re- 
peating the word w has a spectrum determined by the 
stable fixed point value 6(z) corresponding to w. Further- 
more, consider an A^ by A^ matrix H with r^'s given by 
first repeating the word wi (of length Li) Ni times and 
then by repeating the word W2 (of length L2) N2 times. 
As we would expect, in the limit in which A^i, A2, and 
A^ — NiLi + N2L2 all tend to infinity, the spectrum of H 
is given by superposing the spectra of iJ^ (i = 1, 2), where 
Hi is constructed with r^'s given by first repeating the 
word Wi Ni times. This clearly generalizes. In Figs. p|-pl] 

we show the spectrum of the word wi = {+H }, the 

spectrum ofti;2 = {+ + H — }, and the spectrum of the 
word w = wf^W2'^. We see the superposition principle 
at work. 
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FIG. 8: Support of the density of states of all periodic words 
of length 4 or less, over the support of the spectrum of a 
random 4000 x 4000 matrix. 



FIG. 10: Support of the density of states of the periodic word 
{+ + +-}• 








FIG. 9: Support of the density of states of the periodic word 
{+ + —}• 



FIG. 11: Support of the density of states of the composition 
of the two words in Figs. |^ and |l^. 



From this discussion it becomes clear why the spec- 
trum of the matrix H in FZ is so complicated. The se- 
quence {t^i, f2, • ■ ■ Too} is a book written in the binary 
alphabet that, in the mathematical limit N oo, con- 
tains all possible words, sentences, and paragraphs. In 
fact, H contains everything ever written or that will be 
written in the history of the universe, including this par- 
ticular paper. This familiar mind boggling fact accounts 
for the complicated looking spectrum first observed in 
FZ. 

It also explains why numerical studies of the spectrum 
suggest that it is 0— dimensional. Even for N as large 
as 1000 the sequence contains an infinitesimally small 
subset of the set of all possible words, sentences, and 
paragraphs. 



VII. EIGENVALUES ON THE UNIT CIRCLE 

In the ensemble of all books there are particularly sim- 
ple books such that {ri, r2, • • • r^o} consists of a word w 
of length L repeated again and again. In this case, we 
can determine the spectrum explicitly by two different 
methods. 

Let W be the transfer matrix corresponding to w. In 

other words, W — 11 ( f r, 

product is ordered. After repeating the word R times, 
we have 



, where the matrix 



= w 



Ai 
Ao 



(31) 



8 



Diagonalizing W = S 



S, we see 



Ai(z) 
\2{z) 

mediately that A^x is a linear function of A{*(z) and 
A«(z) : 



(32) 



We remind the reader that all quantities in (^2|) are func- 
tions of z. 

The spectrum of H is determined by the zeroes of 
A.rl{z) as i? ^ oo. We note that in this limit the so- 
lution of 



Akl(^) - 



(33) 



does not depend on knowing the detailed form of a and 
p. Indeed, (p^) implies 



or 



(34) 



(35) 



since det W = A1A2 = ±1. In the limit 00, (— /3/q;) « 
tends towards a (z— dependent) complex number of mod- 
ulus unity. Thus, we conclude that 



|A(z)| = l, 



(36) 



namely, that the eigenvalues of W lie on the unit circle. 
This constraint suffices to determine the eigenvalues of 
H. Plugging (p6|), that is A = e'^, into the eigenvalue 
equation 



A^ - (triy)A-t-detT4^ = 0, 



(37) 



we obtain (trT4^) = 2 cos 6* if det = +1 and (trW^) = 
2isin6' if det = —1, which we can combine into the 
single equation 



{tiW) = 2 (det M^) 5 cos 6* 



(38) 



after a trivial phase shift. 

As 9 ranges from to 27r, this traces out the spec- 
trum in the complex plane. As an example, consider the 
4— letter word {+ + H — }, in which case (p8) reduces to 



2z^ = 2icos6'. 



(39) 



This traces out the algebraic curve shown in Fig. [ij, 
which is to be compared to Fig. 

Of course, since H is now translation invariant with 
period L, we can apply Bloch's theorem to determine 
the spectrum of H. Imposing ^k+L = we reduce 

the eigenvalue problem of H to the eigenvalue problem 
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FIG. 12: Support of the density of states of the word {+++ — } 
as given by eq. 

of the L hy L matrix 
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(40) 



One can verify that with a suitable relation between 9 
and (f the eigenvalue equation Det {zIl — h) = becomes 
identical to (p^). 



VIII. CONCLUSION AND OPEN QUESTIONS 

We have seen that the structure of this simple tridiag- 
onal non-Hermitian random matrix possesses an amazing 
richness. This complexity can be understood if one real- 
izes that the spectrum of the random matrix is the sum 
of the spectra of all tridiagonal matrices with a periodic 
subdiagonal obtained by repeating an infinite number of 
times any finite word of length L, weighted by a factor 
1/2^. 

The number of lines does not have the cardinal of the 
continuum. The number of lines is equal to the number 
of words of any length that can be made with a 2-letter 
alphabet; this is a countable number. 

There are many open questions concerning the fine 
structure of the spectrum, such as whether the spectrum 
contains holes in the complex plane (in its domain of def- 
inition). We also have not touched upon the question of 
the nature of the eigenstates. Are they localized or delo- 
calized? Numerical data seems to suggest a localization 
transition. We hope to address these and other questions 
in future work. 
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X. APPENDIX: CYCLIC INVARIANTS 

As remarked in the text, the coefhcients of Ql{z), an 
even polynomial of degree 2L, must be constructed out 
of the cyclic invariants made of wi,'W2,-" ^wl- There 
is presumably a well-developed mathematical theory of 
cyclic invariants, but what we need we can easily deduce 
here. 

For any L, we have the two obvious cyclic invari- 
ants s = X]^=i "-'i ^'^'i P ~ nj'^j^Wj . The number of 
cyclic invariants grows rapidly with L. Apparently dif- 
ferent cyclic invariants can be constructed out of other 
cyclic invariants, for example, ^j) = 

It is easy to work out Ql{z) for low values of L, as 
follows: 

Q2{z)=z^-2sz'' + c£' (41) 
with d = wi — W2, 

Q^{z) = (^3 _ _ 4p ^ ^6 _ 2s_24 _^ ^2^2 _ (43) 

where we have written Q3{z) in a form which shows that 
its roots can be foimd explicitly, 

Qi{z) = z^ - 2sz^ + {s^ + 2k)z'^ - 2sKZ^ + uj^ (43) 



with K = wiws + W2W4 and u> = W1W3 — W2W4, 
Q^{z) = z^°-2sz^+{s^+2k)z^-2skz^+k^z'^~Ap, (44) 

with K = WiWs + Wi'W4 + W2W4 + W2W5 + 

Qq(z) = z^'^ -2sz^° + {s'^ + 2k)z^ -2{sk + p)z^ 

+{k^ + 2sp)z^ -2Kpz'^ + 6"^, (45) 

where 

K = -|- W1W4 + W2W4 + W1W5 + -|- 

+W2W6 + W3We + W4W6 

P = W1W3W5 + W2W4We 

6 = W1W3W5 — W2W4WQ 
and 

Qj{z) = z"-2sz^^ + (s^ + 2k)01°-2(sk + p)z* 

+ 2sp)z^ - 2Kpz^ + p^z^ - Ap, (46) 

with 

K = W1W3 + Wi'W4 + W2W4 + W1W5 + W2W5 + W3W5 

+1^11^6 + W2W^ + WsWq + 'W4W6 + 'W2Wr + WaWj 

+W4Wr + w^w-j 
p = + wiw^wq + W1W4WQ + W2W4We + 

W2W4Wr + W2W5Wr + wzw^w-j. 

The quantities d, k, u, k, p, 6 are manifestly cyclic invari- 
ants. 
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